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ABSTRACT 

Making use of public spectra from Cimatti et al. (2008), we measure for the first time the velocity 
dispersion of spheroid-like massive (M* ~ 10 11 M Q ) galaxies at z ~ 1.6. By comparing with galaxies 
of similar stellar mass at lower redshifts, we find evidence for a mild evolution in velocity dispersion, 
decreasing from ~ 240 kms -1 at z ~ 1.6 down to ~ 180 kms -1 at z ~ 0. Such mild evolution contrasts 
with the strong change in size (a factor of ~ 4) found for these type of objects in the same cosmic 
time, and it is consistent with a progressive larger role, at lower redshift, of the dark matter halo in 
setting the velocity dispersion of these galaxies. We discuss the implications of our results within the 
context of different scenarios proposed for the evolution of these massive objects. 
Subject headings: galaxies: evolution — galaxies: formation — galaxies: structure — galaxies: kine- 
matics and dynamics 
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1. INTRODUCTION 

Recent observations show that the most massive (M* > 
10 11 M Q ) spheroid-like galaxies at z > 1.5, irrespective of 
their star formation activity (Perez-Gonzalez et al. 2008), 
were much smaller (a factor of ~ 4) than their local coun- 
terparts (Daddi et al. 2005; Trujillo et al. 2006, 2007; 
Longhetti et al. 2007; Zirm et al. 2007; Toft et al. 2007; 
Giavalisco et al. 2007; Ravindranath et al. 2008, Cimatti 
et al. 2008, hereafter C08; van Dokkum et al. 2008; 
Buitrago et al. 2008, hereafter B08; Saracco et al. 2009, 
Damjanov et al. 2009, Ferreras et al. 2009). The near 
absence of such systems (r e < 1.5 kpc; M* > 10 n M©) 
in the nearby Universe (< 0.03%; Trujillo et al. 2009) 
implies a strong evolution in the structural properties of 
these massive galaxies as cosmic time evolves. 

Different scenarios have been proposed to explain the 
dramatic size evolution of these galaxies since z ~ 3 
(Khochfar & Silk 2006; Naab et al. 2007; Fan et al. 2008; 
Hopkins et al. 2009). The main difference among them is 
the role of mergers to increase the size of galaxies. Fan et 
al. (2008) support an evolutionary scheme where galax- 
ies grow by the effect of quasar feedback, which removes 
huge amounts of cold gas from the central regions hence 
quenching the star formation. The removal of gas makes 
galaxies to puff up in an scenario similar to the one pro- 
posed to explain the growth of globular clusters (Hills 
1980). In the merging scenario, however, merger rem- 
nants get larger sizes than those of their progenitors by 
transforming the kinetic energy of the colliding systems 
into potential energy. Whereas both scenarios predict a 
strong size evolution for the most massive galaxies, they 
disagree on the expected evolution of the velocity dis- 
persion of the massive galaxies at a given stellar mass. 
The merging scenario basically predicts no evolution (at 
most a 30% since z ~ 3; Hopkins et al. 2009), whereas 
the puffing up scenario predicts central velocity disper- 
sions to be ~ 2 times larger than in present-day massive 
galaxies. 

Constraining the evolution of the velocity dispersion of 
spheroid-like massive galaxies over the last 10 Gyr turns 
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therefore crucial to test the above models of galaxy evo- 
lution, as well as to determine the importance of dark 
matter halos in setting the velocity dispersion of galaxies 
as cosmic time evolves. For these reasons, in this Letter 
we measure for the first time the velocity dispersion of 
such massive galaxies at z ~ 1.6, and compare it with 
a compilation of velocity dispersions for similar galaxies 
at z < 1.2. In what follows, we adopt a cosmology of 
Q m = 0.3, n A = 0.7 and H = 70kms _1 Mpc" 1 . 

2. THE DATA 
2.1. Velocity dispersions at z ~ 1.6 

To determine the typical velocity dispersion of 
spheroid-like massive (M* ~ 10 11 M Q ) galaxies at z ~ 
1.6, we have used the publicly available stacked spec- 
trum of galaxies at 1.4 < z < 2.0 ((z) = 1.63 ± 0.18 
r.m.s. standard deviation) presented in C08 as part of 
the GMASS 1 (Galaxy Mass Assembly ultra-deep Spec- 
troscopic Survey) project. The spectrum consists of an 
averaged spectrum of 13 massive galaxies with a total 
integration time of 480 h. Only 2 of them are classi- 
fied as disks, whereas the remaining 11 galaxies have ei- 
ther a pure elliptical morphology or show a very concen- 
trated and regular shape (see C08). The stacked spec- 
trum is consequently representative of massive spheroid- 
like objects at that redshift. Individual galaxy spectra 
were taken at VLT with FORS2, using the grism G300I 
and a larcsec width slit, providing a spectral resolu- 
tion of FWHM - 13 A in the range AA6000 - 10000 A 
(A/AA ~ 600). Before stacking, each individual galaxy 
spectrum was previously de-redshifted and assigned to 
have the same weig ht in the AA2600 - 3100 A rest-frame 
range. The resulting rest-frame stacked spectrum was 
set to 1 A/pix over the range AA2300 - 3886 A. 

Aimed at computing the velocity dispersion of the 
above stacked spectrum, it is required the use of refer- 
ence template spectra of well known spectral resolution. 
Section l3~2l provides a full description of these data-sets, 
as they are basic input ingredients of the analysis carried 
out in this work. 

1 http://www.arcetri.astro.it/ cimatti/gmass/gmass.html 
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2.2. Velocity dispersions up to z ~ 1.2 

We have compiled from previous work velocity dis- 
persions estimates of spheroid-like massive galaxies at 
z < 1.2. This includes data from van der Wei et 
al. (2005) — updated, when possible, with van der Wei et 
al. (2008) — for galaxies in the ranges 0.6 < z < 0.8 and 
0.9 < z < 1.2, and from di Serego Alighieri et al. (2005) 
for galaxies at 0.9 < z < 1.3. These authors also pro- 
vide stellar masses and effective radii for their galax- 
ies. In both samples the galaxies were classified visually. 
The vast majority of these objects present a prominent 
spheroidal component an are identified as either E or SO. 
When necessary, stellar masses have been transformed 
assuming a Chabrier (2003) IMF. To allow a meaningful 
comparison along the entire redshift interval explored in 
this paper, only galaxies with stellar masses in the range 
0.5 < M* < 2 x 10 n M© have been considered. 

To have a local reference, we have retrieved velocity 
dispersions, effective radii and stellar masses (Chabrier 
IMF) from the SDSS NYU Value-Added Galaxy Cata- 
log (DR6; Blanton et al. 2005; Blanton & Roweis 2007) 
for those galaxies within the above stellar mass range 
having Sersic (1968) indices (in the r-band) n > 2.5. Se- 
lecting objects with n > 2.5 assures that the majority 
of the sources have a Hubble morphological Type T < 
(i.e. ranging from E to SO; Ravindranath et al. 2004). 
We have estimated the average velocity dispersion of the 
above galaxies in two redshift bins: < z < 0.1 and 
0.1 < z < 0.2. We have not tried higher redshifts to 
assure individual signal-to-noise ratios (S/N) larger than 
10. 

3. VELOCITY DISPERSION COMPUTATION 

3.1. The movel program 

To estimate the velocity dispersion of the C08 stacked 
spectrum we use the movel code, available at the RE- 
DUCEME 2 distribution (Cardiel 1999). movel is de- 
signed to characterize the line-of-sight velocity distribu- 
tion of a galaxy spectrum by determining its first two 
moments: the mean radial velocity and the velocity dis- 
persion. Based on Fourier analysis, the code relies on the 
MOVEL and 0PTEMA algorithms of Gonzalez (1993). A full 
explanation of the Fourier techniques and their partic- 
ular implementation in both routines is provided in the 
above work. 

The MOVEL algorithm, an improvement of the classic 
Fourier quotient method by Sargent et al. (1977), is an 
iterative procedure in which a galaxy model is processed 
in parallel to the galaxy spectrum. In this way, a compar- 
ison between the input and recovered broadening func- 
tions for the model allows to correct the galaxy power 
spectrum from any imperfections of the data handling in 
Fourier space. The 0PTEMA algorithm is designed to over- 
come the typical template mismatch problem by comput- 
ing an optimal template for the galaxy spectrum. It is fed 
with a set of representative template spectra which are 
scaled, shifted and broadened according to initial values 
for the mean line strength, the radial velocity, and the 
velocity dispersion. Then, the algorithm finds the linear 
combination of the template spectra that best matches 
the observed galaxy spectrum. This constitutes a first 



composite template which is fed into the MOVEL algo- 
rithm. The output kinematic parameters are then used 
to create an improved composite template and the pro- 
cess is iterated until it converges. This iterative approach 
then provides an optimal template while simultaneously 
computing the radial velocity and velocity dispersion of 
the galaxy spectrum. 

Ideally, the reference template spectra are chosen to 
be representative stellar spectra observed with an iden- 
tical technical setup to that of the galaxy spectra. This 
way, the instrumental broadening is the same in both 
data sets, and the broadening difference between galaxy 
and template spectra, Act, is a direct measurement of 
the galaxy starlight velocity dispersion, er*. Since this 
is not the case for the C08 stacked spectrum, the dif- 
ferent instrumental resolutions of the distinct data sets 
employed in this work must be taken into account. This 
is explained in detail in Sections 13.21 and 13.31 

3.2. The template spectra 

We have restricted the kinematical analysis of the C08 
spectrum to those regions containing strong spectral fea- 
tures from the galaxy starlight. Because of the lack of 
available template spectra of high-enough spectral reso- 
lution covering the full spectral range of interest, we em- 
ploy three different data sets. Two of them allow us to 
analyze the near ultraviolet (NUV) region at A < 3050 A, 
which comprises e.g. several Fen lines at A ~ 2600 A, the 
A2800 Mgn doublet and the A2852 Mgi line. The third 
template set is chosen to provide a complementary anal- 
ysis from the region A > 3250 A, with features such as 
the A3360 NH and A3862 CNO bands. The three sets 
of template spectra, whose properties are summarized in 
Table [TJ are described below: 

i) First, from the CoolCAT 3 database, an HST/STIS 
echelle catalog of normal late-type stars with R ~ 40000 
(Ayres 2005), we have selected eight spectra (from F5V 
to G0III spectral types) covering the range AA2510 — 
3050 A with no gaps in between. They have been re- 
binned at lA/pix to match the linear dispersion of the 
C08 spectrum. Using movel and their corresponding 
high-resolution CoolCAT spectra as reference templates, 
we have derived that the spectral resolution of the re- 
binned templates is ctstis ~ 57±3kms~ 1 at A ~ 2800 A, 
which corresponds to R ~ 2240. 

ii) The second set of template spectra consists of a sub- 
sample of nineteen SSP spectral energy distributions by 
G. Bruzual (private communication). They are based 
on the Bruzual & Chariot (2003) code and the Next 
Generation Spectral Library, NGSL (Gregg et al. 2004; 
Heap &; Lindler 2007), with an average resolution of 
R ~ 1000. The employed SSP models span ages from 
0.3 to 3Gyr with metallicities [Z/H] = -0.4, 0.0, and 
+0.4, hence enclosing the 1 Gyr and solar metallicity es- 
timate of C08 for their stacked spectrum. Again, using 
MOVEL and the above high-resolution CoolCAT spec- 
tra, we have accurately derived the spectral resolution 
of Bruzual's models in the range AA2510-3050A, which 
is f bco3+ngsl ~ 178±8km s -1 , or equivalently R ~ 710. 

iii) Finally, for the red side of the C08 stacked spec- 
trum, we employ a set of eighteen stellar spectra (dwarfs 



2 http://www.ucm.es/info/Astrof/software/reduceme/reduceme.html 3 http://casa.colorado.edu/ ayres/CoolC AT/ 
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TABLE 1 

Velocity dispersion estimates for spheroid-like massive galaxies at z ~ 1.6 



Template spectra Spectral range (A) a A C (A) b cro, tempi (kms 1 ) c cr U: c08 (kms 1 ) d Ac (kms 1 ) ° o"*,C08 (kms 1 ) f 

CoolCAT stars 2510-3050 2780 57 ±3 230 267 ± 20 258 ± 21 

BC03+NGSL SSP models 2510 - 3050 2780 178 ± 8 230 178 ± 23 236 ± 18 

Keck/LRIS stars 3250 - 3880 3565 171 ± 5 179 177 ± 20 236 ± 15 

a Rest-frame spectral range employed for deriving the MOVEL kinematics (Act)* 5 Central wavelength of the rest-frame spectral range c Instrumental 
resolution of the template spectra computed at A c (sec text in Scction l3.2t . d Instrumental resolution of the C08 stacked spectrum computed at A c , 
using FWHM — 13 A and z — 1.6. e Relative velocity dispersion with respect to that of the optimal template spectra/ Velocity dispersion of the 
C08 stacked spectrum of spheroid-like massive galaxies at z~ 1.6, as derived from Equation. 1. 



and giants ranging from GO to K7 spectral types) ob- 
served with LRIS at Keck as radial velocity templates 
for a different program (Beasley et al. 2009). Having 
used the grism 600/4000 and the 1 arcsec longslit, the 
stellar spectra have a FWHM of 4.77 ± 0.14 A (olris ~ 
171 ± 5 kms -1 at 3565 A; R ~ 750) as determined from 
the width of the wavelength calibration arc lines. The 
spectra span the range AA3150 — 5600 A, although only 
the blue side is employed here. 

3.3. Results and robustness of the method 

We have analyzed the stacked spectrum of C08 with 
MOVEL using the three sets of template spectra described 
in Section 13.21 In each case, to constrain the reliability 
of the final results, the procedure has been repeated sev- 
eral times using different initial movel parameters, like 
e.g. different starting velocity dispersions (from 50 up to 
300 kms -1 , in steps of 50 kms -1 ). 

For each single movel solution, the stellar velocity dis- 
persion of the C08 stacked spectrum, ct^coSj is derived 
as 

°*.C08 = V Act 2 + <rg itomp i - ( ^jff y, (1) 

where Act is the broadening difference between the opti- 
mal template spectra and the C08 spectrum (as derived 
from movel), ct tempi and cr^cos are their intrinsic in- 
strumental resolutions, and z is the redshift (~ 1.6). 
Note that, since the individual galaxy spectra of C08 
were previously blue-shifted before stacking, the effec- 
tive instrumental resolution at rest-frame decreases by 
(1 + z). 

For each set of template spectra, Table Q] presents 
Co, tempi, o- 0iC 08, and the mean values of Act and cr*,co8 5 
with their errors accounting for the r.m.s. standard devi- 
ation of the solutions obtained from the different movel 
trials. An error-weighted mean of the three independent 
c*,C08 solutions, which are indeed consistent within the 
uncertainties, gives (cr^cos) = 241 ± lOkms -1 . 

It is worth noting that the above value could be slightly 
overestimated, as redshift uncertainties for the 13 stacked 
galaxy spectra may have enlarged artificially the width 
of the lines. We have also checked that the strong 
A2800 Mgn doublet in the stacked spectrum (~ 8 A, Ta- 
ble 2 of C08) is not significantly broadened by poten- 
tial non-stellar Mgn absorptions in surrounding galactic 
gas clouds. Such absorptions are typically much nar- 
rower (~ 45 kms -1 , as they follow much colder kinemat- 
ics) and ~ 5 times weaker (~ 1.6 A, for typically ~ 14 
clouds; see e.g. Steidel et al. 1994; Churchill et al. 2000). 
Therefore, even if they were present in all galaxies, 



(°*,C08) could be overestimated by only ~ lkms 1 (up 
to ~ 14 kms -1 if all galaxies were Damped-Lya sys- 
tems, which is very unlikely as no Lyman absorptions 
are detected). On the other hand, if the two disk galax- 
ies in the stacked spectrum had lower typical velocity 
dispersions of e.g. 150 kms -1 , the derived (ov.cos) would 
underestimate the true velocity dispersion of early-type 
galaxies by ~ 13 kms -1 . Overall, since the above sys- 
tematics may cancel each other, 241 ± 10kms -1 is prob- 
ably a reasonable estimation of the velocity dispersion of 
M± ~ 10 11 M Q spheroid-like galaxies at z ~ 1.6. 

Figure [T] illustrates representative movel solutions as 
derived from each template set. At the NUV, CoolCAT 
stars (panel a) provide a slightly more broadened solution 
than the NGSL-based SSP models (panel b), probably 
due to the intrinsic differences between the two template 
sets and their effects on the template mismatch problem. 
At longer wavelengths (panel c) , the galaxy spectral fea- 
tures look somewhat blurred, probably by sky line resid- 
uals commonplace of the near-infrared region at observed 
frame. Still, in all cases the solutions reproduce nicely 
the C08 spectrum. 

To further check the reliability of the movel results, we 
carried out Monte Carlo simulations in which we added 
Gaussian noise (to match S/N = 5 to 30 A -1 ) to ar- 
tificially broadened reference template spectra (Act = 
100 — 400 kms -1 , in steps of 50 kms -1 ) from our three 
template sets. For such fake galaxy spectra, we run the 
MOVEL procedure in the same way as we did for the C08 
stacked spectrum. As expected, the accuracy in the final 
solution decreases with the decreasing S/N of the prob- 
lem spectrum. The typical scatter of the derived Act 
solutions — averaged over all the simulations and initial 
MOVEL conditions — are 55, 24, 12, and 8 kms -1 for S/N 
= 5, 10, 20, and 30 A -1 respectively. More importantly, 
for all S/N values, the recovered mean Act is statistically 
consistent with the assumed input value, hence support- 
ing the reliability of both the method and ct^.cos- 

4. DISCUSSION 

The fact that the velocity dispersion of M* ~ 10 n MQ 
spheroid-like galaxies at z ~ 1.6 is ~ 240kms _1 has im- 
portant consequences to understand their evolution: 

• It confirms that high-z spheroid-like massive galax- 
ies are truly massive objects. The unexpected 
strong size evolution of these objects has cast some 
doubts about the reliability of their stellar mass 
estimates, which relies on the assumption that the 
IMF is the same at all redshifts. In fact, an ap- 
propriate change of the IMF with redshift would 
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Fig. 1. — Representative examples of the MOVEL fitting results 
derived for the C08 stacked spectrum and the three sets of template 
spectra: (a) CoolCAT stars, (b) BC03+NGSL SSP models, and (c) 
Keck/LRIS stars. In each case, bottom panels illustrate the COS 
stacked spectrum of galaxies at z ~ 1.6 (red) and the optimal tem- 
plate spectrum (blue, thick line) broadened at the mean Act values 
in Table [T] Top panels show the residuals of the fits. In panel c, 
the region around the O II emission line doublet at A3727 A (dot- 
ted lines) was masked and rejected from the C08 stacked spectrum 
(green) for the fitting computation. 

mitigate the problem of the strong size evolution. 
However, the derived velocity dispersion at z ~ 1.6 
is similar to that of present-day galaxies with M* > 
10 Mq, consequently constraining any potential 
change of the IMF with redshift. 

• It alleviates the problem of understanding how 
massive — and compact — galaxies at high-z can 
evolve through merging since that epoch. An 
extraordinarily high velocity dispersion at high-z 
(e.g. ~ 400 km s -1 ) would have implied that the 



gravitational potential depth of the system would 
be so intense that it would not easily evolve in size. 

To put in context the above result, Figure [2] shows the 
sizes (top panel) and velocity dispersions (bottom panel) 
for the C08 galaxies and the compilation of sphcroid-likc 
galaxies of similar mass described in Section [2] These 
galaxies follow nicely the observed strong size evolution 
found in other independent larger samples where com- 
pleteness effects are accounted (Trujillo et al. 2007; B08; 
dashed line), hence probing that the sizes of the objects 
explored in this paper are typical of average spheroidal- 
like massive galaxies at those redshifts. It is also clear 
that, although massive spheroid-like galaxies have expe- 
rienced a strong size evolution (a factor of ~ 4 since 
z ~ 1.5; Trujillo et al. 2007), the velocity dispersion has 
evolved mildly by a factor of ~ 1.3 in the same redshift 
interval (see Bernardi 2009 for a similar trend at z < 0.3). 
This result is crucial to constrain the different scenarios 
proposed so far to explain the dramatic size evolution 
since z ~ 2, as they disagree in the amount of evolution 
expected in their velocity dispersions. 

In the puffing up scenario of Fan et al. (2008), velocity 

dispersions change as er* cx r c 1 ^ 2 . Using the observed 
size evolution of B08, r c (z) oc (1 + z) -148 , the solid line 
in Fig. [2] illustrates the expected evolution under this 
scenario. It increases with redshift in a way that galaxies 
are expected to double their at z ~ 1.5. The observed 
velocity dispersions are compatible — within the error 
bars — with this scenario up to z ~ 0.7. However, at z > 
1 the discrepancy between theory and data is apparent. 

In the merging scenario of Hopkins et al. (2009), at a 
fixed stellar mass, <t* evolves with redshift as 

a*(z) 1 / r c (0) 

— 7^ K ~7r=i=\/ 7 + — TT' ( 2 > 
cr*(0) Vi+Ty rc ( z ) 

where 7 = (Af na i /i? halo )/(Af*/r c ) is the dark matter 
contribution to the central potential relative to that of 
the baryonic matter at z ~ 0, and i?haio is the effective 
radius of the halo. Again, we have used the B08 size 
evolution as an input into the merging model prediction. 
The grey area in Fig. [2] encloses the expected er* evo- 
lution when 7 varies between 1 and 2. The agreement 
between the merging scheme and the observed evolution 
looks reasonably good at all redshifts. 

Following Hopkins et al. (2009), to explain why a+ 
has only changed weakly since z ~ 2 it is necessary to 
consider that the observed velocity dispersion is driven 
by two components: the baryonic matter and the dark 
matter halo. Both components contribute linearly to 
the central gravitational potential of the galaxy, hence 
cr* oc {M*/r e + M ha io/i?haio)- Assuming that i? ha io 
evolves weakly with time (most simulations show that 
halos build inside-out, so the central potential is set first 
and just the outer halo grows with time), the dark matter 
effect on the central potential of the galaxy (i.e. on a±) 
basically remains unchanged at a fixed M*. However, the 
influence of the baryonic matter on the gravitational po- 
tential has changed strikingly since z ~ 1.5, at present- 
day being ~ 4 times smaller due to the expansion of 
the object. The relative influence of the dark matter on 
setting the inner potential increases with the decreasing 
effect of the baryonic matter. In fact, the data look in 



Mild velocity dispersion evolution of massive galaxies since z ~ 2 



5 



- r 



A von der Wei et al. (2005; 2008) 
• di Serego Alighieri et ol. (2005) - 
■k This work, C08 




1 50 r 



0.0 



0.5 



1 .0 



1 .5 



Fig. 2. — Top Panel: Size evolution of M* ~ 1O 11 M spheroid- 
like galaxies as a function of redshift. Different symbols show the 
median values of the effective radii for the different galaxy sets 
considered in this work (see Section [2}, as indicated in the labels. 
Dashed error bars, if available, show the dispersion of the sam- 
ple, whereas the solid error bars indicate the uncertainty of the 
median value. The dashed line represents the observed evolution 
of sizes r c (z) oc (1 + z)~ 1,48 found in B08 for galaxies of simi- 
lar stellar mass. Bottom Panel: Velocity dispersion evolution of 
the spheroid-like galaxies as a function of redshift, with symbols 
as given above. Assuming the B08 size evolution, the solid line 
represents the prediction from the "puffing up" scenario (Fan et 
al. 2008), whereas the grey area illustrates the velocity dispersion 
evolution within the merger scenario of Hopkins et al. (2009) for 
1 < 7 < 2. See text for details. 



agreement with an almost symmetric influence of dark 
and baryonic matter (i.e. 7 ~ 1) in setting the central 
gravitational potential of present-day objects. 
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